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Abstract 

Using the entire CLEO-c -0(3770) — > DD event sample, corresponding to an integrated luminosity 
of 818 pb _1 and approximately 5.4 x 10 6 DD events, we measure the form factors for the decays 
D° — > p~e + v e and D + — > p°e + v e for the first time and the branching fractions with improved 
precision. A four-dimensional unbinned maximum likelihood fit determines the form factor ratios 
to be: y(0)/Ai(0) = 1.48 ± 0.15 ± 0.05 and A 2 {0)/A l (0) = 0.83 ± 0.11 ± 0.04. Assuming CKM 
unitarity, the known D meson lifetimes and our measured branching fractions we obtain the form 
factor normalizations Ai(0), ^2(0), and V(0). We also present a measurement of the branching 
fraction for D + — > ue + v e with improved precision. 
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The transition rate of charm semileptonic decays depends on the weak quark mixing 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements \V CS \ and \V c d\ [H, an d strong inter- 
action effects binding quarks into hadrons parameterized by form factors. 

In the decays D — » pe + u e , in the limit of negligible lepton mass, the hadronic current 
is described by three dominant form factors: two axial and one vector, A\, A 2 , and V, 
respectively, which are functions of q 2 , the invariant mass of the lepton-neutrino system. 
They are not amenable to unquenched LQCD calculations due to the large total decay 
width of the p meson, but model predictions exist 0, No experimental information on 
these form factors exists. 

The helicity amplitudes for the rare decays B — > V£ + £~ are related at leading order 
in Aqcd/^& to pseudoscalar-to- vector semileptonic transitions [4]. Exploiting one of the 
proposed double-ratio techniques [5j, D — > pe + v e form factors, when combined with those 
of D — > K*e + v e and B — > V£ + £~, can be used to extract \V u b\ from B — > pe + v e . 

The differential decay rate of D — > pe + v e can be expressed in terms of three helicity 
amplitudes (H + (q 2 ), H_(q 2 ), and H (q 2 )) @: 

dT 



dq 2 dcos 9^ dcos 9 e d\ dm 



7Y7V 



J> rGllKdl^^p^Tr^l^W^^p^l + cos^^in 2 ^!^^ 2 ^, ^ 



L 7T7T J 



8(4tt) 4 " cai M 2 D 

+ (1 — cos# e ) 2 sin 2 9 1T \H_(q 2 , m n7r )\ 2 + 4 sin 2 9 e cos 2 9 7T \H (q 2 , m^)] 2 

+4 sin6 l e (l + cos 9 e ) sin 6V cos 6V cosxH + (q 2 , m wn ) H (q 2 , m n7! ) 

—4 sin6 l e (l — cos^e) sin 9 n cos 9 n cosxH-(q 2 , m n7T )H (q 2 , m wn ) 

—2 sin 2 9 e sin 2 9 W cos 2xH + (q 2 , m 7r7r )if_(g 2 , m nn ) , (1) 

where Gp is the Fermi constant, p p is the momentum of the p in the D rest frame, B(p — > tttt) 
is a branching fraction, 9 n is the angle between the ir and the D direction in the p rest 
frame, 9 e is the angle between the e + and the D direction in the e + z/ e rest frame, x is the 
acoplanarity angle between the 7i + n~ and e + v e decay planes, m n7T is the invariant mass of 
the two pions, and BW(m nn ) is the Breit-Wigner function that describes the p line shape. 
Following Ref. 0], we use the relativistic form 



= - 2 „? °ll w„ (2) 



y/m T (p/p ) B(p) 
ml - ml n - zm r(m 7r7r ) B(p ) 

where mo and To are the mass and width of the p meson [§], p is the momentum of 
the pion in the nir rest frame, po is equal to p when m nn = mo, and B(p) is a Blatt- 
Weisskopf form factor given by B(p) = 1/(1 + R 2 p 2 ) 1 ^ , with R = 3 GeV -1 , and r(m 7r7r ) = 
(p/po) 3 (mo/m n7r )T [B(p)/B(p )] 2 . The interference term between a possible s-wave 7T7t com- 
ponent and the p amplitude has not been included in Eq. ([IT) . Its absence is treated as a 
source of systematic uncertainty on the measurement. 
The helicity amplitudes are related to the form factors 

H ± (q 2 ) = MA 1 (q 2 ) T 2^V(q 2 ), (3) 
zm wn -\/q 

~^M*% (4) 
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where Md is the mass of the D meson and M = Mo + m^. Since Ai(q 2 ) is common to all 
three helicity amplitudes, it is natural to define two form factor ratios as 

tv- V(0) andr - M ° ] (5) 

We assume a simple pole form j9[ for A^q 2 ), A 2 ((f), and V(q 2 ), where the pole mass is 
2.01 GeV/c 2 and Mjj*n+\= 2.42 GeV/c 2 |8j for the vector and axial form factors, 
respectively. We have also explored a double-pole parametrization jjjj]. 

We report herein the first measurement of the form factor ratios and absolute form factor 
normalization in D — y pe + u e , and improved branching fraction measurements for these 
decays and D + — y ue + v e . (Throughout this Letter charge-conjugate modes are implied.) 
These decays were studied previously using a smaller CLEO-c data sample [Toj]. The data 
sample used here consists of an integrated luminosity of 818 pb" 1 at the i/j(3770) resonance, 
and includes about 3.0 x 10 6 D°D° and 2.4 x 10 6 D + D~ events. The CLEO-c detector is 
described in detail elsewhere 111!. 



The analysis technique was employed in previous CLEO-c studies [10j,[12|. The presence of 
two D mesons in a DD event allows a tag sample to be defined in which a D is reconstructed 
in a hadronic decay mode. A sub-sample is then formed in which a positron and a set of 
hadrons, as a signature of a semileptonic decay, are required in addition to the tag. The 
semileptonic decay branching fraction Bsl is given by 

n Ntag.SL ftag ^Vtag.SL/e , . 

#SL = — TT = — TT , (6) 

JVtag etag,SL JV tag 

where iV ta g and e tag are the yield and reconstruction efficiency, respectively, for the hadronic 
tag, iVtag,sL and e ta g,sL are those for the combined semileptonic decay and hadronic tag, and 
e = e tag,SL/eta g is the effective signal efficiency. 

Candidate events are selected by reconstructing a D° or D~ tag in the following hadronic 
final states: K + tc~, K + 7i^tt°, and K + 7t~7t~tt + for neutral tags, and K S 7T~ , K + 7t~tt~, 
Kg7r~TT°, K + 7T~7T~7r° , Kg7T~7T~7T + , and K~K + 7i~ for charged tags. Tagged events are se- 
lected based on two variables: AE = E D — E heSbm , the difference between the energy of 
the D~ tag candidate Ed and the beam energy -Ebeam, and the beam-constrained mass 
M bc = (£^ eam /c 4 — |pd| 2 /c 2 ) 1//2 , where is the measured momentum of the D~ candi- 
date. Selection criteria for tracks, 7r°, and Kg candidates used in the reconstruction of tags 



are described in Ref. [13]. If multiple candidates are present in the same tag mode, one 
candidate per tag charge with the smallest \AE\ is chosen. The yield of each tag mode is 
obtained from fits to the M^ c distributions [3]. The data sample comprises 661232±879 
and 481927±810 reconstructed neutral and charged tags, respectively. 

After a tag is identified, we search for an e + and a p~ (V~7r° mode), p° (tt + 7t~ mode), or 
to (7r + 7T~7r° mode) recoiling against the tag following Ref. [13]. A p — > irn candidate satisfies 
\m V v ~ m o| < 150 MeV/c 2 . The combined tag and semileptonic candidates must account 
for all tracks in the event. Semileptonic decays are identified with U = E miss — c|p m i ss |, 
where -Emiss and p m i SS are the missing energy and momentum of the D + meson. If the decay 
products have been correctly identified, U is expected to be zero, since only a neutrino is 
undetected. The resolution in U is improved by constraining the magnitude and direction of 
the D + momentum to be po+ = {E^ C!mi /c 2 — c 2 m|,) 1//2 , and po+ = — p_o- l(j, respectively. 
Due to the finite resolution of the detector, the distribution in U is approximately Gaussian, 
with resolution ~17 MeV for D° — > p~e + v e and D + — y ue + v e and ~8 MeV for D + — y p°e + v e . 
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TABLE I. Signal efficiencies, yields, and branching fractions (Bsl) f° r D° — )> p e + is e , D + — > 
p°e + u e , and D + — > uie + v e , from this work, our previous (prev) measurements [lo| . and two model 
predictions: ISGW2 and FK 0]. All Bsl ar e in units of 10 3 . The uncertainties for e and iVt ag) sl 
are statistical, while the uncertainties for branching fractions are statistical and systematic in that 
order. The efficiencies include the p and u decay branching fractions from the PDG 0]. 



Decay Mode 


6 (%) 


A?tag, SL 




Bsl (prev) 


Bsl(ISGW2) 


%l(FK) 


D° -> p~e+v e 
D+ -> p°e+u e 
D+ —s- u)e + u e 


26.03 ± 0.02 
42.84 ± 0.03 
14.67 ± 0.03 


304.6 ± 20.9 

447.4 ± 24.5 

128.5 ± 12.6 


1.77 ±0.12 ±0.10 
2.17 ±0.12+°;^ 
1.82 ±0.18 ±0.07 


1.94 ± 0.39 ± 0.13 
2.1 ± 0.4 ± 0.1 
1.6±°;5± 0.1 


1.0 
1.3 
1.3 


2.0 
2.5 
2.5 



To remove multiple candidates in each semileptonic mode one combination is chosen per tag 
mode per tag charge, based on the proximity of the invariant masses of the p°, p + , or uj 
candidates to their expected masses. 

The U and invariant mass distributions for D° — > p~e + i> e , D + — > p°e + z/ e , and D + — > 
ue + v e with all tag modes combined are shown in Fig. [TJ The yield for each of the three 
modes is determined from a binned likelihood fit to the U distribution where the signal is 
described by a modified Crystal Ball function with two power-law tails [3] which account 
for initial- and final-state radiation (FSR) and mismeasured tracks. The signal parameters 
are fixed with a GEANT-based Monte Carlo (MC) simulation 15( in fits to the data. The 
background functions are determined by MC simulation that incorporates all available data 
on D meson decays, which we refer to as "generic MC". For D° — > p~e + v e , the backgrounds 
arise mostly from D° — > K*~e + v e , peaking at positive U and modeled with a Gaussian, 
and events with misidentified tags, which are accounted for in the fit by a fourth order 
polynomial. The backgrounds to D + — > p°e + z/ e has its largest contribution from D + — > 
K*°e + i>, K*° — > K~it + , with the peak at higher U due to charged kaons misidentified as 
charged pions, and the peak at lower U from either decay-in-flight kaons or interactions with 
detector material. We categorize the background components according to their shape in 
U and parameterize the overall background shape using combinations of polynomials and 
Gaussian functions. The background shape parameters are fixed in fits to the data, while the 
background normalizations are allowed to float. The signal shapes for the invariant mass 
distributions of the hadronic system are modeled with a Breit-Wigner function, and the 
background shapes are modeled with generic MC. The peaking background for D + — > p°e + z/ e 
arises from D + — > ue + u e , u — > 7r + 7r~. Due to the tag, backgrounds from the non-DD 
processes e + e~ — > qq, where q is a u, d, or s quark, e + e~ — > t + t~ , and e + e~ — > ip(2S)j, are 



negligible |12(] . The signal yields A^ ag SL are given in Table [Q 

The second row of Fig. [1] shows the m^-^o, m^^-, and 777,^+^-^0 distributions with \U\ < 
60 MeV for the three signal modes, respectively. The peaking background at m^-^o ~ 
0.49 GeV/c 2 arises from D° — > K~e + v with K~ — > 7t~tt°. The small background peak at 
m,^- ~ 0.78 GeV/c 2 is due to D + — > ue + u e with u — > tt + tt~ (l6| . 

The absolute branching fractions in Table [U are obtained using Eq. OH]). The signal 
efficiencies e are determined by MC simulation, and have been weighted by the tag yields in 
the data. 

The systematic uncertainties for the branching fractions of D° — > p~e + v e and D + — > 
p°e + z/ e are dominated by uncertainties in the line shape of the p (5.0%), and the non- 
resonant background (—1.5% for D° — > p"e + z/ e and —8.4% for D + — )■ p°e + z/ e ). The un- 
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Invariant Mass (GeV/c 2 ) 

FIG. 1. Fits to the U and hadron invariant mass distributions in data (filled circles with error bars) 
for (a) and (d) D° — > p~~e + v e , p~ — > tt^tt ; (b) and (e) D + — > p°e + v e , p° — > Tr + ir~; and (c) and 
(f) D + — > cje + z/ e , oj —J- 7T + 7r~7r . The solid line represents the fit of the sum of the signal function 
and background function to the data. The dashed line indicates the background contribution. The 
arrows indicate a ±48 MeV region around the Kg mass, which has been removed for display. 



certainty due to the line shape of the p is estimated by (1) requiring \U\ < 60 MeV and 
fitting the distribution, (2) varying the selection criterion \m nn — m \ < 150 MeV. The 
uncertainty due to the non-resonant background is obtained by performing a form factor 
fit, with an additional interfering non-resonant D — > 7T7re + z/ e (s-wave) component modeled 
following Ref . [17| , then integrating over the kinematic variables to recalculate the branching 
fractions. The unknown form factors in D + — > ue + v e are the dominant uncertainty in its 
branching fraction (3.0%). The remaining systematic uncertainties include the track and 7T° 
finding efficiencies, positron and charged hadron identification, the number of tags, the no- 
additional-track requirement, the shape of the signal and background functions, and the MC 
FSR and form factor modeling. These estimates are added in quadrature to obtain the total 
systematic uncertainties on the branching fractions: _~io.o%; 4.1%, for D° — > p~e + v e , 

D + — > p°e + i> e , and D + — > ue + u e , respectively. 

A form factor analysis is performed for D — > pe + u e . We calculate the energy and mo- 
mentum of the neutrino using E v = E miss and |p„| = -E m i SS , because i? m i SS is better measured 
than |pmi SS |. Without ambiguity, the four kinematic variables (g 2 , cos^, cos 6^ \) are 
measured with resolutions of (0.021 GeV 2 /c 4 , 0.020, 0.048, 0.024) for D° ->- p~e + v e , and 
(0.013 GeV 2 /c 4 , 0.013, 0.037, 0.019) for D + -> p°e + u e . 

A four-dimensional maximum likelihood fit in a manner similar to Ref. 18( is performed 
in the space of q 2 , cos9 n , cos 6^, and \- The technique makes possible a multidimensional 
fit to variables modified by experimental acceptance and resolution taking into account cor- 
relations among the variables. The signal probability density function for the likelihood 
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FIG. 2. Projections of the combined p~ and p° data (points with statistical error bars) and the fit 
(solid histogram) onto q 2 , cos# e , cos 9^, and \- The dashed lines show the sum of the background 
distributions. 



function is estimated at each data point using signal MC events by sampling the MC distri- 
bution at the reconstructed level in a search volume around the data point, then weighting 
by the ratio of the decay distribution for the trial values of ry and r 2 to that of the gener- 
ated distribution. The search volumes are one tenth the full kinematic range of each of the 
four dimensions. Large MC samples are generated to ensure that each search volume has 
sufficient statistics. The background probability density function is modeled using events 
from the generic MC. Due to the low statistics of the background in the generic MC, we 
reduce the four dimensional space to lower dimensional subspaces. Due to the correlation 
between q 2 and cos^e, the two subspaces are chosen to be (q 2 , cos9 e ) and (cos^, x)- The 
background normalization is fixed in the fits to the values measured in the determination of 
the branching fractions. 

Using the above method, a simultaneous fit is made to the isospin-conjugate modes 
D° ->■ p-e + u e and D+ ->■ p°e + u e . We find ry = 1.48 ± 0.15 and r 2 = 0.83 ± 0.11, with a 
correlation coefficient pv2 — —0.18. The confidence level of the fit is determined to be 5.0% 
by comparing the negative log-likelihood from the data to the distribution from toy MC fits. 
Fig. |2] shows the q 2 , cos# e , cos^, and x projections for the combined p~ and p° data and 
the fit. We also make fits to the two modes separately. The results are consistent. We note 
that the difference between the data and the fit projection for cos 8 n might be due to s-wave 
interference. 
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We have considered the following sources of systematic uncertainty in the form factor 
measurement. Our estimate of their magnitude are given in parentheses for ry and r 2 , 
respectively. The uncertainty associated with background modeling (0.01, 0.02) is estimated 
by changing the normalization of the three largest background components by a factor of 
two in each semileptonic mode. The uncertainty due to imperfect knowledge of the p line 
shape (0.01, 0.02) is estimated by modifying the p line shape by increasingand decreasing 
the population of signal MC events below and above the nominal p mass |8j by 20%. The 
uncertainty due to non-resonant background (0.01, 0.02) is obtained by repeating the fit with 



an additional interfering non-resonant D — > -K-Ke + v e component (s-wave) following Ref. 17 
The procedure for extracting the form factor parameters is tested using the generic MC 
sample, from which events are drawn randomly to form mock data samples, each equivalent 
in size to the data sample. When backgrounds are absent, the measured form factor ratios 
are consistent with the input values. In the presence of background, a small statistically 
significant shift is observed. Its magnitude is taken as the uncertainty due to possible bias 
in the form factor fitter (0.03, 0.02). The uncertainty associated with the unknown q 2 
dependence of the form factors (0.03, 0.02) is estimated by introducing a second pole [jjj]. 

Adding all sources of systematic uncertainty in quadrature, the final result is ry = 1.48 ± 
0.15 ± 0.05 and r 2 = 0.83 ± 0.11 ± 0.04. Using \V cd \ = 0.2252 ± 0.0007 obtained using 
CKM unitarity constraints |8j and the lifetimes t d o = (410.1 ±1.5) x 10 _15 s and t d + = 
(1040 ± 7) x 10 _15 s jsj, we combine our form factor ratio and branching fraction results to 
obtain Ai(0) = 0.56 ± 0.0llg;g§, A 2 (0) = 0.47 ± 0.06 ± 0.04, and V{0) = 0.84 ± 0.09±^. 



Our branching fraction results are compared to previous measurements [10J, with which 
they are consistent, and theoretical predictions in Table HI The results are consistent with 
isospin invariance: 2T^D^ p v e +u) = 0-09io!o2- I sos P m symmetry is not expected to 

be exact due to p° — uj interference [l6j. Theoretical predictions from the ISGW2 model [2| 
and a model (FK) which combines heavy-quark symmetry and properties of the chiral La- 
grangian [jjj], are also listed in Table HI The branching fractions for ISGW2 are obtained 
by combining the partial rates in Ref. 0] with \V c d\ and td from PDG 0. Our branching 
fraction results are more consistent with the FK predictions than ISGW2. 

The FK model predicts A^O) = 0.61, A 2 (0) = 0.31, and 1/(0) = 1.05. These values are 
compatible with our form factor measurements. No other experimental form factor results 
on these decays exist. Our values of ry and r 2 are very similar to the current PDG average 
of D + -»■ K*°e + i> form factor ratios r v = 1.62 ± 0.08 and r 2 = 0.83 ± 0.05 



In summary, we have made the first measurement of the form factor ratios and absolute 
form factor normalization in D — >■ pe + u e , and improved branching fraction measurements 
for these decays and D + — > ue + u e . Our branching fractions are consistent with our previous 
measurements but with improved precision. The form factor measurement in D — > pe + u e is 
the first in a semileptonic Cabibbo-suppressed pseudoscalar-to-vector transition. 
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